Physical Review D 78, 025006-7 (2008) 




DMA [iMECC] UNICAMP 





FERMION-FERMION BOUND STATE CONDITION FOR 
g : SCALAR EXCHANGES 

^ ■ 

^SJ . Stefano Dc Leo and Pictro RotcUi 

Department of Applied Mathematics, University of Campinas 



00 



CD 



2 



PO Box 6065, SP 13083-970, Campinas, Brazil 
deleo @ime. unicamp. br 

Department of Physics, University of Lecce and INFN Lecce 

PO BOX 193, CAP 73100, Lecce, Italy 

rotelli@le.infn.it 

Submitted April, 2008 



Abstract. The condition for the existence of a bound state between two fermions ex- 
changing massive scalars is derived. For low scalar mass, we reproduce the scalar field 
model result. The high scalar mass result exhibits a somewhat different inequality con- 
dition. 
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g : I. INTRODUCTION 



Since the experimental discovery of the mass of the neutrinos [1] a legitimate question has been posed. 
Is there a possibility of a bound state between weakly interacting particles such as an electron and a 
neutrino? If the particles involved where nonrelativistic the answer would be immediate and negative. 
From numerical studies of the Schrodinger equation [2] the existence of a bound state produced by a 
Yukawa (attractive) coupling 



OO 

O 

> 

g^exph/.r] 
4:77 r 

has been found to be 

^>0.84-^, (2) 

where ^ is the exchanged particle mass and m is the reduced mass. A related derivation is the use 
of a surrogate to the Yukawa potential, the Hulthcn potential, Vh, which approximates the Yukawa 
potential for small r, 

^"^''^""4^exp[2Air]-l " 

Our choice of is made so that the terms r^^ and r° in a scries development about r = are 
identical to the Yukawa potential. The Schrodinger equation with the Hulthcn potential can be solved 
analytically [3] and the existence of a bound state yields a similar result to that above, i.e. 

f^>^. (4) 
47r m 

High mass exchanges would necessitate extremely strong couplings, obviously unphysical for weak 
interactions where ^/m > 10" [4]. However, if one considers the relativistic corrections to the 
Schrodinger equation one encounters the well known Darwin term [5] 

^vV(r), (5) 
8 m 
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which for a Yukawa potential V^{r) yields 

1 



ti^V^{T)+AnglJ{r)] . (6) 

O lit' 

The first term above can be summed with the potential contribution to yield an overall amplification 
factor 

(7) 

8 nv- J 

This is what has been called Yukawa coupling amplification [6]. Here the effect must be small to 
comply with the very nature of correction terms. However, if one where so bold as to assume this 
amplification for high ^/m one would invert the resonance condition i.e. 

8 m 

^ > , 8 

which allows bound states even for the weak interactions. The problem, theoretically, now shifts to 
determining the resonance condition for high mass exchanges in a more rigorous manner. A method 
has been introduced and applied in field theory [7, 8] . It consists of confronting the lowest ladder 
contributions (box and crossed) to the scattering amplitude at rest, with the tree diagram contribution 
(also in the rest frame). The requirement that the sum of the former be equal or greater than the tree 
contribution reproduces exactly the Hulthen condition for low ///to within a scalar-scalar model with 
scalar particle exchanges. In this paper it will also be shown to be also valid in the case of fermion- 
fermion (f-f) interacting with scalar exchanges. More recently [9], the scalar model calculation was 
extended to the high /i/m limit (in either limit, approximations or numerical calculations are needed). 
The high /i/m result was even more restrictive than the Hulthen inequality (|4]), i.e. it required even 
larger gl^^, specifically 

^>^/ (luJ^+^-±J^Xup] , (9) 

with p = rrii/m^. However, it was noted that since the Klein-Gordon equation lacks a Darwin term 
correction there is no reason to expect Yukawa amplification. In this paper, we essentially repeat our 
low and high p/m limits for f-f interacting via scalar exchange. This case does contain a Darwin term 
identical to that of the well known electrostatic case although some additional corrections also exist. 

In the next Section, we illustrate the model and reduce the first order ladder contributions to 
a single integral in d|fc| = Ak. In Section HI, we perform the small p/m limit and reproduce the 
Hulthen inequality (|4]). In Section IV, we perform the high p./m limit. We propose a phenomenological 
expression for the k integral based upon numerical simulations. In Section V, we draw our conclusions. 

II. THE FERMIONIC MODEL 

In the center of mass system and for forward scattering (see Fig. I), the Feynman rules [4] for the 
amplitudes of the box (□) and crossed box (x) diagram yield 

. ^ . . 2 /• ^ ^l'^^ (-P) { [E, (p) + E,{p)]^,^^+m,} u'f'^ i-p) 4"^ jp) i^ + m,) uj"'^ (p) 
^P^-'^^^^J {2ny D^{p)D,{p)Dl{p) 

(10) 

d^fc u^^\-p) {^ + m, + [E,{p) - E,{p)] l„}v!{'\-p) ui'\p) {f: + m,)ui''\p) 



M^{p)=iglgl 
with 



(27r)* Dnp)D,{p)Dl{p) 



Xs 



u['Kq) = ^E,,,{q)+m,^, [ ^.g ^ | {s ^ 1,2) , Xi^[q} , ^'=[l 

Ei,2{q) + rui 
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The denominators factors are, 

D?{p) = Elik) - [fc„ - E,{p) - E,{p)] - ie , 
ip) = eI (fc) - [ fco + E, ip) - E,{p)] - le , 
D,{p) ^ Elik) ^ kl - le , 

D„{p) ^ Elik - p) - [fco - E,ip)] ~ ie , (11) 

where 

At threshold (p - 0), 

d*fc (kg + TO2)(2mi + TO2 " fco) 



M''''{0) = i {2g^g^^mi'm.,f 5rr' 5., 



{2^Y D^{0)DMDm 

d*fc (fco + 2 TOi — ma) (fco + JTla) 



(12) 



The poles in the lower half complex fco plane are at 

fc°i = Ei{k) + nil + , k^^^ ^ Ei{k) - + , fco, 2 = E^ik) , fco,o = Eo{k) + . 

The box and crossed box diagrams give the following fourth-order contribution to the invariant scat- 
tering amplitude 



A1°(0)+ A^X(O) = i^M^^^Srr'S.,J dU- g [i??(fc)+i?r(fc)" 

^ (^g^gW^^ J Srr' S,.,, J^"° dk e [i?°(fc)+i?X(fc)] . (13) 

Below by Es we intend Es{k) and by W and A we intend TOi -|- ma and ?7i2 — TTii respectively. A 
simple calculation shows that the explicit formulas for the residues in the fco-plane for the box and 
the crossed box diagram are respectively 

i?°(fc) = [2m,{m,-E,)-k-'] / {4W E, {E, + m,) [i^l' - 2m, {E, + m,)]"} , 

R°{k) = [fc" - 2mi(£;2 + m2)] / { 4 1^ i;2 (£^2 - "^2) [^f" + 2 7712 (S2 - 7^2) ]'} , 

i?„°(fc) = [2(£o + 2m2)(2mi-So)] [{E,, - m,) B C + (E, + m,) AnC - AnB] / [A^^B^C] 

+ 2 {m,-m,-E,) / [AaBC^] , (14) 

with Aq ~ 2 Eo rrii — fi'^ , B ~ —2 E„ ?7i2 — and C = 2 Eg, and 

i?f = [e + 2m.,{E,+m,)] / {A E, {E, - m,)[^l'' + 2 m, {E, ~ m.,)f] , 
= - [fc" + 2mi (£2+^2)] / {4A£;2 (^2 - m2) [m' + 2m2 (i;2 - 7712)]'} , 
=2{Eo + 2m,){Eo + 2m,)[{Eo + m,) B C + {Eo + m,) A^C^A^B]/ [AlB'C] 

+ 2 {E, + W) I {A^BC^\ , (15) 

with = — {2E„m, +/i^). It is to be noted, and can be used in calculation, that the residues of 
the box and crossed residues are related by 

^i'^2,o = - -R°2,o[™i ^ -rrii] . (16) 

However, care must be used when applying this symmetry because for example ^/ml + m, = 2 m,, 
while, under m, — !■ — m,, {—m,Y ~ ^n, = 0^—2 m,. The rule of thumb is that square root factors 
should be left as such before applying such symmetries. 

Before passing to the actual calculation of the small and large ^/m results, we must discuss two 
important technical questions. The first is the question of the convergence of the fc integrals. The 
second is the feature of real pole contributions in some of these residue integrals. 
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• Convergence. 

Individually, the leading residues terms yield divergent integrals, both linear and logarithmic. This 
was not the case for the scalar model [8]. However, when summed, the divergences cancel, specifically 
in the limit k oo, 

16 ml mlk-^R^ = -'4 + ' "^^"^^ " ^ "^^'^^ 



W Wm,k 
□ ^ m'l [jj.^ + 2 17111712 — 3 ml)ml 




16 TO? ml k'RV = + — - 

' ' ' W Wm^k 

16 ml ml fc- = + A + ^^'^ + ' "^^^^^^"^^ + "^p'^ ' ^ + O f 1 

TOi TOa K \ fc^ 

Consequently, 

16TO>^Fi?° = -^^ + of-^^ and 16m>?Fi?>< =+^$^ + f-^ 

Both these results lead to convergent integrals, however, when summed, the leading terms again cancel 
and finally 

16TO>^F(i?°+i?x) = -^5^ + o(l) , (17) 

which is a highly convergent integrand. Notice that this leading order result is symmetric under 
TOi <-> 7712. We have not specified which mass, TOj or ma, is the lower mass and the Feynman diagrams 
are clearly symmetric under the interchange TOj TO2. Any final results must therefore be symmetric 
under this symmetry. This feature may be used as a test of all of the following results. 



• Poles. 

By explicit observation the quadratic term in the denominator of i?P vanishes at /i^ = 2toi {Ei +mi). 
Poles also occur in the expression for when Aq — 0, i.e. at fx' = 2toi Eq. Both of these conditions 
correspond to the same value of fc, which we indicate by kp, 

No other residues have poles. Thus, R^ and R\^ exhibit double and single poles on the real axis 
at kp. However, when summed all pole contributions cancel. This is demonstrated in some detail in 
the Appendix. The cancellation of the double pole is simple to show. That of the single pole which 
receives a contribution from R^ and four contributions from R^ , one from each term in the last line 
of Eq. (|14p . is more cumbersome to sec. However, it must be proved since it would otherwise dominate 
the large /i/to calculation, and radically change our conclusions. 



III. THE EXCHANGE OF SMALL MASS SCALARS 

For incoming fcrmions with mass TOi and TO2 interacting by the exchange of a scalar with mass 
fj. <C TOi,2, and R^ contribute to the invariant scattering amplitude only for value of fc <IC TOi_2 
(indeed of the order of /i). In this small /i limit, wc may use the approximation 

Ei,2 = J + ml 2 « TOi,2 + TT" — ■ 

We note, as an aside that for small fi (^ TOi^2) there are no poles on the real axis. Now it is easy to 
show that 

R?/R? = 0[(M/m)«] « 1 . 
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Whence in the rest of this Section i?^ will be neglected. The other residue contributions yield 

e Rl 
k^Rl 
F Rl 



>□ ^ 


2 miTOj 


1 


+ 


1 fc^ 


mi fc*' 


-•2 ^ 


W 






2W^ E^ 


m, 


X ^ 


2 miTOa 


1 




1 fc^ 


mi fc"* 




A 


K 




2A 


m, A ' 




2 m-ifno 


1 


+ 


1 fc^ 


mj fc* 




^+ A 






2A ^ 


^ mi A i?„^ ' 


>□ ^ 


3 fc" 




A 




/i^A P 


'0 ^ 


' + 4 £f " 


2mim.2 Eg 


2 mi ma i?o 


X ^ 


3 fc" 




W 






'0 




2niim2 Eg 


2 mi ma Eg 



Thus, 



i?°(fc)+i?''(fc) 



1 

2W 



1 

771, 



Et 



1 



1 

m. 



and, by making use of the elementary integrals 

Fdfc 16 fi 



TT 



4^ 

TT 



wc find that 



' dfc 



Et 



37r 



2m 



fcM/s _ 16 

TT m 
" 2 ^ 



Fdfc 



mi 



(19) 



= 1 



57r 1 



Comparing now this fourth-order total scattering amplitude, 



^2 2 2 



TT w ^J.^ 

with the one boson exchange amplitude (tree diagram) 

9i92 



1 



8 iriim^ 



4 mima 



(20) 



(21) 



(22) 



we find that the fourth-order amplitude is greater or comparable to the second-order amplitude when 



1.92 ^ }J_ 



47r 



777-1 ma 



(23) 



which, to the leading order, reproduces exactly the Hulthcn inequality, where gl^^ ~ 9i92- Wc have 
explicitly calculated and exhibited the correction term in the above inequality, and wc will refer to 
this factor in our conclusions. 



IV. THE EXCHANGE OF HIGH MASS SCALARS 

The high ^/m limit is more difficult to treat and we rely upon numerical tests of the following 
expressions. We have three masses in our calculation of A^^'^ so if we consider an adimensional 
expression, it can only be a function of /i/mi and /i/ma or alternatively of 



and 



mi 
m. 



Indeed, 



dk k" 



R°ik) + R''{k) =F{uj,p) , 



(24) 
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and this can be tested numerically. Now we try to parameterize Mf^'^ by a form derived in the scalar 
model case. We write, 



dk fc' 



a 



In I 



1 



l-p' 



hip 



(25) 



The value a = \ reproduces the scalar model result. This phenomenological form has been tested for 
a wide but limited range of uj and p values, specifically for 



w = 10% 10' , 10" and p = 2 , 10 , 50 



In the following Table 





p 


Phcn/Nimr 


10" 


2 


.995 


10" 


10 


.989 


10" 


50 


.978 


10' 


2 


1.005 


10' 


10 


1.000 


10' 


50 


.993 


10" 


2 


1.012 


10" 


10 


1.008 


10" 


50 


1.003 



we give the comparison of phcnomcnological/numerical (Phcn/Num) results for a best fit value of a, 

a = 0.663 . (26) 

We see that to within a few per cent the agreement is good. Wc could of course improve the comparison 
if we included a constant term In (3 in the brackets which could correspond to e renormalization of 
the logarithmic terms. However, we consider this an excessive fincss. The important point is that the 
large /x/m behavior is similar to the scalar model result. The high p, resonance inequality thus reads 



27r 



a TOim2 



1 



(27) 



V. CONCLUSIONS 

We have applied in this paper a field theoretic approach to the determination of the coupling strengths 
needed for the existence of a fcrmion-fermion bound state via scalar boson exchanges. For low p/m, 
we again find the Hulthcn inequality [3,8] as seen in the scalar model. For high p/m, we obtain an 
even more restrictive condition (P7)) . a result again similar to the scalar field model [9]. The similarity 
between the scalar field model and this calculation suggests that the bound state inequality condition 
depends essentially upon the exchanged particles rather than the incoming ones. This was by no 
means obvious since the numerators of the residues are different in the two cases. Indeed at first sight 
the fermion-fermion model seemed to yield divergent results as a simple power count of the fc-integral 
suggests. We have shown in this paper that the individual divergence contributions cancel. We have 
also shown that the real pole contributions to A^'-' also cancel both for the double and single poles. 
Again it is not clear if this would happen with say vector particle exchanges and it must be said that 
a contribution from a simple pole would completely alter our high p/m results. For the existence of 
a relativistic bound state such a contribution could even be desirable. 

There is however a problem with our results for small p/m and the arguments based upon the 
relativistic corrections to the Schrodinger equation mentioned in the introduction. The Dirac equa- 
tion with a scalar potential contains a Darwin term as does the better known electrostatic case [5]. 
This lead us to expect, at least for small p/m (nonrclativistic) a coupling amplification. We have 
purposefully kept the 0{p'^/m^) corrections in the small p/m case and as can be seen in the result 
([23|) the corrections terms correspond to a coupling deamplification. The coupling constants must be 
somewhat increased to compensate the correction terms. This result is consistent with the tougher 
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large ///m inequality. We predict that the Huhhen inequahty is a lower limit inequality for any /i/jTi. Is 
this disagreement between our field theory calculation and the nonrelativistic reduced mass equation 
serious? This may well be a matter of opinion but some observations are in order: 

- The Hulthen inequality is not exactly in agreement with the Yukawa numeric inequality. So, we 
have a formal discrepancy even neglecting the relativistic correction terms; 

- The higher order Feynman diagrams cannot be parameterized by a simple Yukawa potential. 
However, the Coulomb potential works admirably well for Hydrogen like atoms except for one of the 
supreme successes of field theory, the Lamb shift. Unfortunately, we known of no direct way to derive 
the potential bound state spectrum from field theory; 

- It must also be remembered that not all the fourth order Feynman diagrams have been calculated. 

Nevertheless, we remain troubled by this result. At the very least, we must moderate any expectations 
for a weak interaction calculation in which intermediate vector particles are exchanged. We expect 
the same low /i/m inequality (except perhaps for the correction term) but hope for a very different 
high /x/m result. 

Our results have one physical consequence, we predict that weak interacting fermion-fermion (or 
scalar-scalar) particles cannot produce a bound state simply by Higgs boson exchanges [1]. It is our 
intention to tackle the full weak interaction case in the near future. 
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APPENDIX: THE POLE CONTRIBUTIONS 

In this Appendix, we calculate the pole contributions of i?^ and at 



For convenience, we define the functions F{k), G{k) and H{k) by 

k^R^{k) = F{k) , (28) 

3 

= G(fc) + ^H„(A:) , (29) 

where 

G{k) = 2k'{E„ + 2m,)(2m,-E,){E,^m,) / (A^B C=) , 
H,{k) = 2k^E„ + 2m,){2m,-E,){E, + m,) / {AaB^G') , 
H,{k) = -2e{Ea + 2m,){2m,-E,)/ {AaB C) , 
Hsik) = 2P (mi -£^o) / {AaB C^) . 

The first pole terms in the MacLaurin series of these functions are 

( — ■ 2) — 2 

where F (kp) is the coefficient of {k — kp) in F{k) and so forth. 
Now for the double pole, we find that the only two contributions are 

(kp) = ~d-\kp)^ (2K- m;)(4^.^. + M-) (30) 

whence their sum cancels. 

The single pole contributions can be written as 

{f'-(M,G'-'(M.«'-(M}^3j^ {/;-.»;-.'•;"■'} . m 

and in the Table we list the factors in graph brackets above as a series in even powers of fi. For 
example, 

4 2 

= -2W fi' - {m, + A)W ^ ~ m,W ^ . 
^ 2 nil 8 nil 

The important point is contained in the last line of the Table. All single pole contributions also cancel. 
Thus, in conclusion, there are no real axis poles in fc' (i?° + i?^). 
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/i* / 2 mi 




/i" / 3 mlm2 


(-1) 

Jp 


-2VF 









(' 1) 
9p 


m^ + 2W 


2 Aiy - 2 toJ + 3 TO1TO2 


2 7712 A 


-(7772 + 


hpA 




2 — 3 rriimz 


7712(2 TTli — A) 




, (-1) 





-m,W 


-AW^ 




n-p-s 





2m,W 


AW 





Sum 















Table. The coefficients of powers of /i^ in the factors fp \ ' and hip * for the single 
contributions. These factors are defined in Eq. ([?T|) . 
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Fig. 1 

[E^{p),-p] [E,(p)+E^(p)-ko,-k] [E,(p),-p] 
> • ^ • ^ 

i/ [ko- E^{p),k-p] A 

^ • > • ^ 

[E^{p),+p] [ko,k] [E^{p),+p] 

box diagram 



[E,{p),-p] [ko + E,(p)-E^(p),k] [E,(p),-p] 
^ ^ ^» ^ 




^ 9^ ^ • ^ 

[E,ip),+p] [ko,k] [E,ip),+p] 

crossed box diagram 



Fig. 1. The fourth order box and crossed box diagrams in a ferminonic field model evaluated in the center of 
mass frame for scattering in the forward direction. 



